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We present the gravitational waveforms computed in ab initio two-dimensional core collapse supernova mod¬ 
els evolved with the CHIMERA code for progenitor masses between 12 and 25 M©. All models employ multi¬ 
frequency neutrino transport in the ray-by-ray approximation, state-of-the-art weak interaction physics, rela¬ 
tivistic transport corrections such as the gravitational redshift of neutrinos, two-dimensional hydrodynamics 
with the commensurate relativistic corrections, Newtonian self-gravity with a general relativistic monopole cor¬ 
rection, and the Lattimer-Swesty equation of state with 220 MeV compressibility, and begin with the most 
recent Woosley-Heger nonrotating progenitors in this mass range. All of our models exhibit robust explosions. 

Therefore, our waveforms capture all stages of supernova development: 1) a relatively short and weak prompt 
signal, 2) a quiescent stage, 3) a strong signal due to convection and SASI activity, 4) termination of active 
accretion onto the proto-neutron star, and 5) a slowly increasing tail that reaches a saturation value. Fourier 
decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO for 
Galactic events across the range of progenitors considered here. The fundamental limitation of these models is 
in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models, which 
are underway. 


I. INTRODUCTION 

Core collapse supernovae (CCSNe) are recognized as the 
most energetic explosions in the modem Universe. As a re¬ 
sult of the collapse and the subsequent explosion of a massive 
star (M > 9Mq), gravitational binding energy of the pro¬ 
genitor (100 B = 10®^ ergs) is released in the form of neu¬ 
trino radiation ( 99%) and the kinetic energy of ejected mate¬ 
rial observed in the electromagnetic spectrum (~ 1%). Post¬ 
collapse asymmetries in fluid motion along with anisotropic 
emission of neutrinos also generate a strong burst of gravita¬ 
tional waves (GWs). The presence of three different types of 
radiation associated with CCSNe makes them ideal objects for 
multi-messenger astronomy [1]. Simultaneous observations 
of the multi-messengers from CCSNe can reveal not only the 
details of the supernova mechanism [2] but also may shed 
light on fundamental properties of neutrinos [3]. Gravitational 
radiation signals are particularly interesting among them be¬ 
cause they are the only signals that provide deep insight into 
the multidimensional dynamics of the supernova core. The 
uniqueness of the information carried by GWs was recognized 
quite early (for comprehensive reviews of GWs from CCSNe, 
see [4-6]). 
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However, CCSNe are extremely complex and physically di¬ 
verse phenomena that involve an intricate interplay of gen¬ 
eral relativistic gravity, hydrodynamics, and neutrino trans¬ 
port, and thermonuclear kinetics, on short time scales. The 
complexity of the problem requires state-of-the-art numerical 
simulations for a quantitative analysis of the processes taking 
place in the supernova core. This is a great challenge in and of 
itself. It is currently known that realistic simulations of CC¬ 
SNe are multidimensional and include hydrodynamics, self¬ 
gravity, and neutrino transport with the complete set of neu¬ 
trino weak interactions, preferably in full general relativity but 
with at least approximate general relativistic corrections to all 
three. Recent simulations of other groups that have provided 
gravitational waveforms met some of these requirements [7- 
10 ]. 

The first simulations of GW emission from core collapse 
supernovae were limited to rapidly rotating stellar core col¬ 
lapse [11-14]. These simulations did not require a large 
amount of computational resources in order to produce the 
strongest part of the GW signal at bounce and, therefore, their 
primary waveforms are of relatively short duration (^SO¬ 
SO ms). More recent predictions of the signal from rapidly 
rotating collapsing cores [9, 10, 15-24] use two- or three- 
dimensional simulations, conformally-flat or full general rel¬ 
ativity, along with, in some cases, a deleptonization prescrip¬ 
tion for the stellar core, developed to reproduce the results of 
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a full neutrino transport treatment during the collapse phase 
[25] and a neutrino leakage scheme for the post-bounce evo¬ 
lution [22]. For the relatively short simulation time associated 
with these models, this approach is quite adequate. However, 
most supernova cores likely do not rotate rapidly [26] to pro¬ 
duce a strong GW signal. 

Another class of CCSN simulations has focused on the GW 
signatures of hydrodynamic instabilities in the post-bounce 
phases - in particular, convection inside the proto-neutron star 
(PNS), neutrino-driven convection in the post-shock region 
[27-31], and the standing accretion shock instability (SASI) 
[32-38]. These models require that we simulate at least sev¬ 
eral hundred milliseconds of post-bounce time during which 
time all instabilities develop and achieve their nonlinear state, 
and potentially saturate, and during which time explosion also 
develops. In this case, a far more realistic treatment of neu¬ 
trino transport is crucial. Although there has been recent 
signihcant progress (e.g., see [39]), multidimensional CCSN 
models with realistic neutrino transport and full general rel¬ 
ativity are not yet available. This is in part due to the shear 
computational cost associated with such models that cannot 
be paid on present-day supercomputing platforms. There are 
a few methods widely used to approximate neutrino transport 
in two- or three-dimensional models that do not require as 
much computational power: a parameterized approximation 
for neutrino heating and cooling [4, 7, 18, 31], gray neutrino 
transport schemes [8, 40], the IDSA method [41], and finally, 
but most notably, multi-group neutrino transport (e.g., multi¬ 
group flux-limited diffusion) in the “ray-by-ray-plus” approx¬ 
imation [42]. Previous studies [7, 8, 43, 44] using one of 
these methods have established the fundamental structure of 
the waveform, which consists of three major components: a 
prompt-convection signal that lasts about 70 ms; the strongest 
part of the signal, due to neutrino-driven convection and the 
SASI (~200 ms); and a signal due to the revived shock’s ex¬ 
pansion, which produces a growing offset in GW amplitude. 
Uniqueness of the “ray-by-ray” method is that it allows us to 
estimate the contribution of anisotropic neutrino emission to 
the total GW signal [8, 43^5]. 

In this paper, we present the gravitational waveforms from 
four non-rotating axisymmetric (two-dimensional) relativis¬ 
tic models evolved with the neutrino-hydrodynamics code 
Chimera beyond one second of post-bounce time and ex¬ 
plosion. We compare our results with both GW predictions 
from our previous simulations (A-series) [44] and the results 
of other groups [8, 9]. We discuss all of the features of the 
GW signals from the four models initiated from four different 
progenitors with masses 12, 15, 20, and 25 Mq. All of our 
waveforms are available for download from ChimeraSN.org. 

Our paper is organized as follows. A brief description of 
the code and the model setup is given in Section 2. In Section 
3, we describe the methods used for extraction and analysis of 
the GW signals. In Section 4, we provide gravitational wave¬ 
forms for all of our models and the results of our analysis. In 
Section 5, we summarize our investigation and draw conclu¬ 
sions. 


II. CODE DESCRIPTION AND MODEL SETUP 

We analyze the GW emission in four two-dimensional sim¬ 
ulations performed with the neutrino-hydrodynamics code 
Chimera [46] . Chimera consists of five major mod¬ 
ules: hydrodynamics, neutrino transport, self-gravity, a nu¬ 
clear equation of state, and a nuclear reaction network. 

Chimera solves the equations of Newtonian hydrody¬ 
namics but takes into account some effects of strong-held 
gravity by means of an “effective potential” [47, 48]. The 
gravitational held is computed by multipole expansion [49]. 
We replace the Newtonian monopole component with a GR 
monopole [48]. It has been shown by Muller et al. [50] that 
this approximation works very well for slow rotation. The 
neutrino transport module solves the energy-dependent neu¬ 
trino moment equations for all neutrino havors using an up¬ 
dated version of multi-group Hux-limited diffusion (in the ray- 
by-ray-plus approach of [42]) with a Hux limiter that has been 
tuned to reproduce Boltzmann transport results [51]. The ray- 
by-ray-plus method is able to produce angular variations in 
the neutrino radiation held, which, in turn, generates low fre¬ 
quency GW signals (Figure 1). 

In total, we evolve four different non-rotating, non- 
perturbed, axisymmetric models (designated B12-WH07, 
B15- WH07, B20-WH07, and B25-WH07, corresponding to 
zero-age main sequence progenitors of 12, 15, 20, and 25 Mq 
[ 52]) on a spherical-polar mesh consisting of 512 non-equally- 
spaced, adaptive radial zones and 256 uniformly-spaced angu¬ 
lar zones. In radius, the grid covers 3x 10"^, 2 x 10“^, 2.1 x 10^, 
and 2.3 x 10"* km, respectively. In angle, the grid goes from 0 
to TT. 

All presented models were simulated using the Lattimer- 
Swesty [53] equation of state (EoS) with a bulk incompress¬ 
ibility modulus K = 220 MeV for p > lO^^gcm”^, which 
is capable of supporting the maximum observed neutron star 
masses of ^2 Mq [54, 55], and an enhanced version of the 
Cooperstein EoS [56] for p < lO^^gcm”^ where nuclear sta¬ 
tistical equilibrium (NSE) applies. 

All four progenitors are evolved beyond 1 second after core 
bounce. All models exhibit shock revival and the development 
of neutrino-driven explosions. The evolution beyond 1 sec¬ 
ond captures all important phases of the supernova dynamics 
that pronounce in the GW signals. Three of the models de¬ 
veloped clear prolate shock morphologies, while the 20 Mq 
model develops an approximately spherical, off-center shock 
as the explosion begins, and then becomes moderately prolate 
at ~600 ms after bounce. The morphologies of the explo¬ 
sions are reflected in the gravitational waveforms. The explo¬ 
sion geometry of the 20 Mq model has reduced not only the 
model’s explosion energy relative to the 15 and 25 Mq models 
but also the total energy emitted in the form of GWs. 

III. GRAVITATIONAL WAVE EXTRACTION 

We consider the slow-motion, weak-fleld approximation for 
post-processing extraction of the GW signal from our simula¬ 
tions. 
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A. GWs Produced by a Time-Dependent, Mass-Quadrupole 
Moment 

We consider only the lowest order terms in the retarded ex¬ 
pansion of the mass-quadrupole formula [57]. The transverse- 
traceless part (TT) of the gravitational strain is given by 



where D is the distance from the source to the observer and 
the mass quadrupole l 2 m {t — is defined as 

hm = / rooY2*mr^dV. ( 2 ) 

Here rpo is the (OO)-component of the linearized stress-energy 
tensor [58]. In the weak-field limit. Too is approximated by 
the rest-mass density of matter. The tensor spherical harmon¬ 
ics <!>) (d and (p are the angular coordinates of the ob¬ 

server’s frame (O-frame)) are defined in the Appendix. The 
amplitude A 2 m in the gravitational strain is 


Optimal design of gravitational wave detectors requires some 
knowledge of the expected waveforms, the corresponding 
frequency spectra, and the total energy emitted by possible 
sources [59]. Thus, any extraction method should decrease 
the numerical noise as much as possible. Most numerical dif¬ 
ferentiation methods amplify the numerical noise built into the 
simulation data. To avoid this, A 2 m is usually computed by 
reducing the order of time derivatives of l 2 m.' 

dN2m , dl2m ... 

Mm = , where N 2 m = —j—- (4) 

at at 

Following [60, 61] and [62], the quadrupole signal can be ex¬ 
pressed in terms of a volume integral depending only on the 
density, velocity, and the gradient of the gravitational poten¬ 
tial; 


where r, t) and ip are the spherical coordinates in the source 
frame (S'-frame), is the component of velocity in the same 
frame. 

Further reduction of the time derivative is conventionally 
done by using the momentum equation [60, 61]. However, 
this has to be carried out carefully. The Euler equation in¬ 
cludes stress terms. Therefore, in order to replace dt{pv) with 
balancing stress terms, we need to take into consideration all 
possible contributions to the stress terms, such as pressure, 
gravity, anisotropic neutrino forces acting on the fluid, effec¬ 
tive viscosity that may present in the finite differencing of the 
momentum equation, etc. To avoid this issue, we have decided 
to compute N 2 m for all time steps and to numerically evaluate 
dN 2 m/dt to obtain A 2 m- Numerical algorithms for comput¬ 
ing first-order derivatives introduce far less numerical noise 
than those for higher-order derivatives. In order to determine 
the detectability of the GWs, we calculate the characteristic 
GW strain for a given frequency / [63] using 
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where dEgw if) /df is the GW energy spectrum and A 20 if) 
denotes the Fourier transform of ^20 it)' 


00 

A2oif)= I A2oit)e-^^^f^dt. (8) 

— 00 


The stochastic nature of GW signals from CCSNe prompts 
the use of short-time Fourier transform (STFT) techniques to 
determine the frequency of a signal as it changes over time 
[7]: 


00 

STFT{A2o(f)}(T,/) = J A2o(f)iZ(f-T)e-*2-/‘df (9) 

— 00 

where Hit — t) is the Hann window function [64]. 


= J (5) 

Using the continuity equation in the integrand in Eq. (5), one 
can replace the density time derivative. Integrating by parts 
and omitting the surface contribution, we find the resulting 
integrand 
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B. Gravitational Waves Produced by Anisotropic Neutrino 
Emission 

Besides aspherical mass motion, any other sources with 
non-zero quadrupole moments will produce GW emission. 
One of these sources is the anisotropic radiation of neutrinos 
from the hot PNS. The theoretical derivation of the GW sig¬ 
nal produced by a distant anisotropic point source of neutrinos 
was first published by [65]. Muller end Janka [31] were the 
first authors to implement this formalism. Kotake et al. [4] im¬ 
proved the formalism and made it more suitable for numerical 
evaluation of GW signals. 

The transverse-traceless part of the gravitational strain, 
hj^, from neutrinos is given by [31] 
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Here H is the solid angle in the 5'-frame, and the vector 
is the direction of neutrino emission whose components are 
given with respect to the O-frame. The tensor {riinjY"^ is the 
transverse-traceless part of the second-rank symmetric tensor 
Uirij with respect to the observer’s z-axis (the z-direction is 
defined as the direction connecting the source and the ob¬ 
server). The angles 9 and (j) define the direction of neutrino 
emission with respect to the O-frame. The other factor in the 
integrand, dL/dfl, is the “direction-dependent neutrino lumi¬ 
nosity” given in the S'-frame. 

In the case of axisymmetry, both /i+ and hx components 
of the gravitational wave signal vanish for an observer on the 
symmetry axis, and the GW signal with the maximum ampli¬ 
tude will be detected by an observer in the equatorial plane 
[31]. For such an observer, the gravitational strain is given by 
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studies of CCSNe based on models that explode [7, 8, 43, 44]. 
The GW signal passes through four distinct phases: 1) A 
prompt signal: an initial and relatively weak signal that starts 
at bounce and ends at about 80 ms post-bounce. 2) A qui¬ 
escent stage that immediately follows the prompt signal and 
ends somewhere about 125 ms after bounce. 3) A strong sig¬ 
nal, which follows the quiescent stage and is the most ener¬ 
getic part of the GW signal. This stage ends somewhere be¬ 
tween 350 ms and 400 ms after bounce. 4) A “tail,” which 
starts just before the end of the strong phase, at about 300 
ms after bounce and consists of a low-frequency component 
with increasing amplitude. This “tail” tends to rise during the 
simulations, but not monotonically. The signal produced by 
anisotropic neutrino emission exhibits a low-frequency “tail”- 
like behavior during the entire simulation. 

The frequency evolution in Figure 1 shows how the charac¬ 
teristic frequency of the signal changes during the simulation. 
The initial frequency ~500 Hz of the early signal drops to 
100-200 Hz (the quiescent stage). The strong phase of the 
signal is characterized by a steady increase of the peak fre¬ 
quency, which reaches its maximum value of ~700-800 Hz 
at 500-600 ms. Then, the frequency slowly decreases during 
the “tail” phase. 


B. Early Signal 


T t a! I\ /I • a/ cos^ r?'- sin^ 0'sin^ 

T'(r?,Vj) = (l + sin 0 cos^j)- ^ . 2 . 2 , ■ 

cos^ w + sin i)' sin (p' 

( 12 ) 

Since our models are axisymmetric, we can simplify 
41 {'O'T ip'). Integrating over ip' [4]: 


41 {d') = sin0' —TT + / dp' 
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IV. WAVEFORM ANALYSIS 

Gravitational waveforms covering all phases of GW emis¬ 
sion and based on non-parameterized supernova explosions 
were first reported in the work of [44]. Here we present wave¬ 
forms obtained in a new series of supernova simulations per¬ 
formed by our group [66], and a comparison with our previous 
results. 


A. Qualitative Description of the GW Signals 

GW signals from both matter and neutrino sources for 
the models B12-WH07, B15-WH07, B20-WH07, and B25- 
WH07 are provided in the left panels of Figure 1 , and the evo¬ 
lution of the signal frequency is given in the right panels of the 
same figure. Qualitatively, our gravitational waveforms have 
all of the key features of waveforms described in previous 


The early GW signal is produced by entropy-driven con¬ 
vection in the PNS during the first 10 ms after bounce. Unlike 
in our previously published models (the “A-series”) [44], in 
the models presented here (the “B-series”)[66] we do not ob¬ 
serve the high-amplitude, low-frequency contribution to the 
prompt signal due to the deflection of infalling matter through 
the shock. As a result, the peak amplitudes of the prompt sig¬ 
nals in all the B-models decrease by a factor of 2—4 relative to 
their values in the corresponding A-series models (see the in¬ 
sets of Figure 2). The results presented here are in agreement 
with the conclusions drawn by [8], who estimated the con¬ 
tribution at the shock in their semi-analytic approach under 
the assumption that the non-radial component of the veloc¬ 
ity is negligible. For our A-series runs, this assumption was 
not valid. In that series, we set the lateral velocities of the 
fluid to zero above the shock during a short time after bounce 
for numerical reasons, which led to a coherent deflection of 
collapsing matter when passing through shock, with a sudden 
and commensurate change in the 0-velocities of the fluid ele¬ 
ments. One can see in Figure 2 the effect is more pronounced 
for a lesser progenitor mass. This happens because the rela¬ 
tive contribution of the lateral velocity to the total velocity of 
the fluid is higher for less-massive progenitors (and, therefore, 
with lower radial velocities). In the B-series models, asym¬ 
metries in the angular velocities above the shock are obtained 
given commensurate asymmetries in the gravitational poten¬ 
tial. This, in turn, reduces the jump in the angular velocities 
across the shock and excludes the presence of the coherent de¬ 
flections seen in the A-series models and described above, as 
illustrated in Figure 2. 

As was correctly pointed out in [8], prompt convection can- 
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FIG. 1. Left panel: The GW signals generated by both matter (solid gray) and neutrinos (dashed black) for our 12-25 Mq models. Right 
panel: Spectrogram showing a normalized value of /df as a function of frequency and time after bounce. The blue line tracks the 

evolution of the peak frequency of the signal. 
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post-bounce time (ms) 


FIG. 2. Comparison of the waveforms obtained in our A- and B- 
series CHIMERA simulations for the 12, 15 and 25 Mq models. 


not be the only source of the quasi-periodic signal in this 
emission phase, due to the difference in time scales of prompt 
convection (~30 ms) and the prompt signal (~80 ms). After 
prompt convection has ceased, we do observe a diminishing 
but nonzero signal amplitude. The reason is that the shock 
expands quickly during the first 80 ms after bounce. There¬ 
fore, the perturbed matter is accumulated in the fast-growing 
volume behind the shock. Although convection has become 
less pronounced after 20 ms, integration over an increasing 
volume in part compensates for this. 

Further clarification of the contributions to the prompt GW 
signal would require simulations with a finer grid resolution 
that would be capable of capturing all hydrodynamic instabil¬ 
ities (i.e. prompt convection, turbulent cascades, acoustic and 
vorticity waves, etc.) that are presented in this phase. 


C. Strong Signal Phase 

The beginning of the strong signal phase coincides with the 
onset of SASI activity. It has been shown in previous stud¬ 
ies [7, 43, 44] that the strong signal is actually produced by 
the combined effect of SASTinduced downflows impinging 
on the PNS surface and the subsequent deceleration of the 
matter at the PNS surface and the convection inside the PNS. 
The low-frequency component arises from the modulations in 
the shock radius as the SASI develops and evolves. The high- 
frequency component is generated when the SASTinduced ac¬ 
cretion flows strike the PNS surface (Figure 3). It is clear from 
the analysis of the contributions to the strain from r < 50 km 
and r > 50 km that the PNS convection, deceleration of the 
accreting matter at the PNS surface, and neutrino-driven con¬ 
vection in the gain region contribute significantly. 




FIG. 3. Top: The entropy distribution for the B15-WH07 model 
inside the PNS at 228 ms after bounce. Down flows onto and con¬ 
vective activity inside the high-density region produce the strongest 
GW signal. 

Bottom: The GW waveforms, Dh+ vs. time, showing the contribu¬ 
tions of three regions: r < 50 km, r > 50 km and r > 500 km. The 
latter region shows the contribution due to the shock expansion. 


The stalling of the shock manifests itself in a relatively short 
quiescent phase, which lasts about 20-30 ms and which is 
followed by a strong signal produced by the development of 
neutrino-driven convection and the SASI. 


The shock modulations affect the kinetic energy of the ac¬ 
cretion flows and, consequently, the amplitude of the GWs 
generated when these flows hit the PNS surface. The signal 
structure during the strong signal phase in both B12-WH07 
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and B15-WH07 is similar to that in the corresponding A- 
series models. However, this is not the case for B25-WH07 
and A25-WH07. The beginning of the strong signal phase in 
A25-WH07 is ~50 ms behind that in B25-WH07, which in¬ 
dicates an earlier development of neutrino-driven convection 
and SASl activity in the latter model. The peak amplitude in 
B25-WH07 is twice as large as it is in A25-WH07. 

The peak frequency of the signal grows almost linearly 
from 100 Hz up to 1000 Hz during the strong signal phase 
(right panels of Figure 1). We see the same trend in frequency 
evolution, with a similar slope, in the M15 model from Muller 
et al. [8], which is the closest to our B15-WH07 model. 


D. Explosion Phase 

All of our GW signals end with a slowly increasing tail, 
which reflects the (linear) gravitational memory associated 
with accelerations at the prolate outgoing shock. The notice¬ 
able decrease of the high-frequency component of the ampli¬ 
tude during the explosion phase (most pronounced for model 
B12-WH07 at 520 ms) is due to the cessation of active accre¬ 
tion onto the PNS surface (Figure 4). The time of the ces¬ 
sation coincides, within a width of the STFT window, with 
the time when the frequency reaches its maximum value, for 
all of our models except B20-WH07 (Figure 1). B20-WH07 
has a different explosion morphology. A single downstream 
is formed in all of our models except B20-WH07 in the early 
SASI phase. This downstream produces the local large am¬ 
plitude spikes in the GW strain by its deceleration at the PNS 
surface. The downflow also induces the 1=2 mode of the mass 
distribution deep in the PNS, which enables high-frequency 
PNS convection to contribute to the GW signal. Thus, PNS 
convection is responsible for the high-frequency component 
of the GW waveform. Termination of the single accretion 
stream leads to a significant decrease in both the frequency 
and the amplitude of the GW signal. In B20-WH07, mul¬ 
tiple downstreams are formed during the SASl phase. This 
prevents the establishment of a more precise correlation be¬ 
tween the changes of the accretion flow and the associated 
changes of the waveform amplitude and peak frequency. The 
typical frequency in B20-WH07 starts to decrease when the 
first accretion downflow detaches from the surface of the PNS 
(~500 ms) while other downstreams continue to perturb the 
PNS and thus support the high-frequency and the amplitude of 
the B20-WH07 signal (Figure 5), until the moment when the 
last accretion down flow becomes detached from the PNS sur¬ 
face (^630 ms). After the cessation of accretion, the GW sig¬ 
nal in all of our models is essentially generated by the shock 
only. The tails continue to rise until they reach their satura¬ 
tion values at 700-1000 ms, depending on the model and its 
prolateness. 

The total emitted GW energy is shown in Figure 6. The 
values of the GW energy emitted in the B-series models pre¬ 
sented here are very close to what we predicted in the A-series 
models presented in [44]. Due to the “anomalous”evolution 
of model B20-WH07, we do not observe a simple correlation 
between the progenitor mass and the total energy emitted in 




post-bounce time (ms) 


FIG. 4. The entropy distributions for the B12-WH07 model are 
shown at two post-bounce times: just before the accretion down¬ 
stream finally detaches from the surface of the PNS (top panel) and 
just after the detachment (middle panel). This detachment produces 
a clear signature in the gravitational waveform. 

Bottom panel: The amplitude of the high-frequency component of 
the GW signal sharply decreases at the time (~530 ms) when the 
active accretion rate onto the PNS surface ( ~50 km) drops signifi¬ 
cantly (right scale). Further variations of the amplitude of the signal 
are mainly the result of the interaction between the outgoing shock 
and the infalling matter. The blue lines locate the significant change 
in the rate of change of the accretion rate, at ~450 ms. 


gravitational waves. That is, the GW energy emitted is not 
solely dependent on the progenitor mass. It is also a func¬ 
tion of the explosion dynamics - in particular, the number and 


Accretion Rate (M/s) 
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B20-WH07 



FIG. 5. Entropy distribution snapshots of the B20-WH07 model 
that show multiple down streams observed during the early accretion 
phase (top panel)', a single accretion down stream at a significantly 
later post-bounce time (middle panel)', the moment of final detach¬ 
ment of the final accretion down flow from the surface of the PNS 
(bottom panel). This explains why the time when the peak frequency 
of the signal reaches its maximum (~500 ms) doesn’t correlate well 
with the time when active accretion onto the PNS finally ceases. In 
B20-WH07, we begin with multiple accretion streams. The peak 
frequency begins to decline as the number of accretion streams is re¬ 
duced, not as the single accretion stream detaches from the PNS, as in 
the other models. Note also: The significant change in the compact¬ 
ness of the PNS due to neutrino radiation emphasizes the importance 
of neutrino transport and general relativity in supernova simulations 
and to the accurate prediction of the associated gravitational wave¬ 
forms. 


characteristics of the accretion streams that form during the 
pre-explosion and explosion phases. The GW energy emitted 
does increase monotonically with progenitor mass for mod¬ 
els B12-WH07, B15-WH07, and B25-WH07, which share the 
same explosion morphology, but the results are very different 
for B20-WH07 given the difference in its explosion morphol¬ 
ogy and the resultant difference in the evolution of its accre¬ 


tion streams. 


T-'-1-'-1-'-r 




FIG. 6. Top: The gravitational waveforms for all of the models pre¬ 
sented here. Bottom: The energy Eq'n radiated in the form of gravi¬ 
tational waves as a function of time. 


As shown in Figure 6, almost all of the GW energy is emit¬ 
ted between 200 ms and 700 ms after bounce for all of our 
models, and we do not observe significant contributions to 
i?Gw in any model after 700 ms, when the low-frequency com¬ 
ponent of the GW signal becomes dominant. The jumps in 
the emitted gravitational wave energy correspond to abrupt in¬ 
creases in the accretion rate onto the PNS surface. This is eas¬ 
ily seen by correlating these jumps with the counterpart spikes 
in the gravitational waveforms in Figure 6 that are the direct 
result of such accretion. Our predictions for the gravitational 
wave energy emitted are 2-3 times higher than those based on 
the general relativistic models G11.2 and G15 of Muller et al. 
[8]. This is largely due to the fact that the high-amplitude, 
GW tail contributes significantly to the GW energy emitted. 
In turn, such large-amplitude tails are not produced in models 
that do not explode robustly. 
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E. Signals Produced by Neutrino Emission 

The GW signals from neutrino emission in the B-series 
models are similar to the signals in their A-series counterparts, 
with some variation associated with the nonlinear stochastic 
nature of multidimensional supernova models. The ampli¬ 
tudes of all of the GW signals from neutrino emission are 
slightly negative just after bounce, until ~ 180-220 ms, de¬ 
pending on the model, and then increase dramatically, becom¬ 
ing positive throughout the rest of the simulation. This is the 
result of the formation of a stable accretion down flow in the 
central region of the grid (60° < 0 < 120°) due to an active 
interplay between neutrino-driven convection and the SASI at 
this time. The cold dense matter in the formed down flow ab¬ 
sorbs neutrinos more efficiently than the matter in the polar 
regions. As a consequence, the neutrino luminosity is more 
intense in the polar regions, which makes the amplitude of 
the GW signal positive (Eq. 13). The situation is more in¬ 
tricate in B20-WH07, due to the presence of multiple down 
streams, but the general trend is similar to the other models. 
Note that the amplitude of the neutrino-generated GW signal 
is much larger than the amplitude of the matter-generated GW 
signal; however, neutrino-induced GW signals have relatively 
low frequencies (/ < 20 Hz) for the canonical setup of cur¬ 
rent gravitational wave detectors. Nevertheless these signals 
may be detectable using the “non-traditional” approach pre¬ 
sented in [67, 68]: three-mass experiments may allow one to 
measure a permanent displacement of the test masses due to 
the (linear) memory effect of gravitational waves. 


F. Detectability of the Signals 

Figure 7 compares the GW strain spectra h^harif) of 
our models with the broadband design noise levels of 
advanced-generation GW interferometers, assuming a source 
distance of 10 kpc. Most of the detectable emission is 
within ~ 100-800 Hz, with the level increasing from (^2 to 
10)xl0“^^ Hz“^/^. A Galactic event (at 10 kpc) appears 
to be well detectable by the upcoming generation of detec¬ 
tors. The peaks at ~550-750 Hz are due to a cumulative 
effect of high-frequency convection inside the PNS and the 
deceleration of downflows at the PNS surface. In general, the 
peak frequencies of all of our models presented here are lower 
than those seen in the A-series models because of the low- 
frequency contribution of the late signal ( >600 ms). Though 
all of the peaks lie in a relatively narrow frequency interval, 
one can see that the peak frequency tends to decrease with 
increasing progenitor mass. 

Our GW predictions for the B15-WH07 model can be com¬ 
pared to the M15 model of Muller et al. [8] given that both 
groups implement similar treatments of the neutrino transport 
and GR corrections to the gravitational field, and include es¬ 
sentially the same physics in their models. The two groups 
are in agreement with regard to the time scales of the different 
(pre-explosion) GW phases, and the amplitude of the prompt 
and strong GW signals. They differ, however, in their predic¬ 
tions for the peak in the GW spectrum, which is at ~1000 Hz 


in the M15 model versus 650 Hz in model B15-WH07. This 
difference likely arises due to the presence of a strong explo¬ 
sion in our simulation that considerably decreases active ac¬ 
cretion onto the PNS after 400 ms post bounce and, conse¬ 
quently, the peak frequency of the GW signal (Figure 7). 


-B12-WH07(xl0^ 
-BI5-WH07(xl0^ 
- B20-WH07 
-B25-WH07(xl0) 


'h 


IE-20 


IE-21 


IE-22 


IE-23 ' 


IE-24 



Frequency (Hz) 


FIG. 7. Top: The GW energy spectra for all of our B-series mod¬ 
els. In order to better differentiate the curves, the spectra have been 
rescaled by a factor of 10“^, 10“^, and 10 for models B12-WH07, 
B15-WH07, and B25-WH07, respectively. 

Bottom: The characteristic GW amplitudes, hchar, for all of our B- 
series models, plotted against the approximate noise thresholds for 
Advanced LIGO, Advanced Virgo, and KAGRA, at a source loca¬ 
tion of 10 kpc. 

Dashed, colored lines mark the peak frequencies of the correspond¬ 
ing models on both spectrographs. 


V. SUMMARY AND CONCLUSIONS 

Based on four ab initio axisymmetric explosion models 
[46] for non-rotating progenitors with masses 12, 15, 20, and 
25 Mq, we studied the GW signals of core collapse super¬ 
novae from the early post-bounce phase through the fully- 
developed explosion phase. Unlike in our earlier studies [44], 
which were truncated at ^500 ms after bounce, here we pro¬ 
vide the complete (up to 1 s) gravitational waveforms for all 
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four models. This is the hrst time complete signals have been 
computed in the context of ab initio explosion models. This 
is particularly important for the signatures in frequency space, 
whose accurate determination requires a full temporal evolu¬ 
tion. 

Our models qualitatively confirm the four-phase picture of 
GW emission seen in our previous studies, and by others us¬ 
ing parameterized models, with an early quasi-periodic sig¬ 
nal, a quiescent phase of several tens of milliseconds, a strong 
stochastic GW signal lasting until some fraction of a second 
after the onset of explosion, and a low-frequency tail. 

Given that we tracked the full dynamical and GW evolution 
beyond 1 s after bounce, we were able to follow the transitions 
between the four stages of the GW signal and, in particular, 
between the last two stages, at the moment of cessation of ac¬ 
tive accretion onto the PNS surface. This is clearly seen in the 
gravitational waveforms, both in the behavior of the strain as a 
function of time and in the evolution of the peak frequency of 
the GW signal, and is especially manifest in the B12-WH07 
model. The evolution of the peak frequency of the signal, 
which declines monotonically after accretion has stopped, is 
clearly evident in all four cases. Moreover, the peak frequency 
in the GW energy spectrum and the characteristic strain is in¬ 
versely related to progenitor mass. On the other hand, the 
total GW energy emitted does not exhibit a simple correlation 
with progenitor mass. For models B12-WH07, B15-WH07, 
and B25-WH07, which all exhibit the same prolate explosion 
morphology, the GW energy emitted increases monotonically 
with progenitor mass. However, model B20-WH07, with its 
more spherical explosion, emits signihcantly less energy in 
GWs than models B25-WH07 and B15-WH07, which indi¬ 
cates that the energy emitted in GWs is a gross function of 
two parameters: progenitor mass and explosion morphology. 

Our findings clearly demonstrate that a complete prediction 
of gravitational waveforms and spectra is possible only for 
a ab initio simulation with a fully developed explosion. Both 
the peak frequency of the GW signal and the GW energy emit¬ 
ted depend on the duration of the GW signal. In the latter case, 
the GW energy emitted does not saturate until 400-700 ms af¬ 
ter bounce, depending on the progenitor. 


In addition to carrying out our simulations past 1 s after 
bounce, differences between the signals provided here and 
those published in our earlier work [44] were obtained for the 
early, hrst-phase signal and, as discussed here, were the result 
of different model specifications - in particular, constraining 
the pre-shock flow to be spherically symmetric or not. The 
results described here were obtained in models that do not im¬ 
pose such a constraint, which is artificial. 

A major shortcoming of the work presented here is our 
imposition of axisymmetry. Complete GW waveform pre¬ 
dictions, including both the /i+ and polarizations, in the 
context of state-of-the-art three-dimensional supernova simu¬ 
lations are required. The expected differences in the stellar 
core hydrodynamics in three dimensions case [9, 24, 69-78] 
will have an impact on the predicted waveforms. In particular, 
changes in the geometry of the accretion down flows and the 
behavior of turbulence may influence the amplitude and en¬ 
ergy of the GW emission. We have already reported a success¬ 
ful development of supernova explosion in the 3D CHIMERA 
simulation [79]. Within a month, we will present an overview 
of the simulation [80] followed by the paper presented the 3D 
gravitational waveforms. 

In addition to the leap to three dimensions, other “dimen¬ 
sions” of the problem - e.g., the use of different nuclear 
equations of state, progenitors (especially non-spherical 
progenitors (see [81]), etc. - need to be considered, as well, 
in the context of late-time three-dimensional models. 
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ions and conclusions expressed herein are those of the authors 
and do not represent the National Science Foundation. The 
authors would like to acknowledge fruitful discussions with 
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Appendix A: Tensor Spherical Harmonics fim 


In this appendix, we write explicitly the tensor spherical 
harmonics. We follow the notation of [4]. 


1. Explicit Form of Wim and Xim 


2. Tensor Spherical Harmonics 


Let us define the tensor spherical harmonics (in the O- 
frame) using the Wim and Xim functions 


In general, the Wim and Xim functions are defined by 
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Xim{9,cl)) = 2^ - cot 9 '^ Yim{9,cl)), (A2) 


where Yim is the usual scalar spherical harmonics of order I 
and degree m. 

For I = 2 (quadrupole approximation), we have the follow¬ 
ing explicit expressions: 
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Here the first row (column) corresponds to 9 (</)), and a is a 
normalization factor. Notice, this tensor is trace-free, and the 
diagonal components correspond to the -f mode, while the 
off-diagonal components correspond to the x mode. 

The normalization is fixed by the following relation: 


J dn iflm)AB ^ = Sll'Smm', 

(A14) 

where A,B,C,D = 9,(1), and the metric on the 2-sphere, 
is 
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The normalization factor for Z = 2 is 
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